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ABSTRACT 

We study the abundance of silicon in the intergalactic medium by analyzing the statistics of Si IV, 
C IV, and H I pixel optical depths in a sample of 19 high-quality quasar absorption spectra, which 
we compare with realistic spectra drawn from a hydrodynamical simulation. Simulations with a 
constant and uniform Si/C ratio, a C distribution as derived in Paper II of this series, and a UV 
background (UVB) model from Haardt & Madau reproduce the observed trends in the ratio of Si IV 
and C IV optical depths, Tgiiv /tctv • The ratio tsuy/tcw depends strongly on Tcrv, but it is nearly 
independent of redshift for fixed tctv, and is inconsistent with a sharp change in the hardness of the 
UVB at z « 3. Scaling the simulated optical depth ratios gives a measurement of the global Si/C ratio 
(using our fiducial UVB, which includes both galaxy and quasar contributions) of [Si/C] = 0.77 ±0.05, 
with a possible systematic error of ~ 0.1 dex. The inferred [Si/C] depends on the shape of the UVB 
(harder backgrounds leading to higher [Si/C]), ranging from [Si/C]~ 1.5 for a quasar-only UVB, to 
[Si/C]~ 0.25 for a UVB including both galaxies and artificial softening; this provides the dominant 
uncertainty in the overall [Si/C]. Examination of the full tsuy/tcw distribution yields no evidence 
for inhomogeneity in [Si/C] and constrains the width of a lognormal probability distribution in [Si/C] 
to be much smaller than that of [C/H]; this implies a common origin for Si and C. Since the inferred 
[Si/C] depends on the UVB shape, this also suggests that inhomogeneities in the hardness of the UVB 
are small. There is no evidence for evolution in [Si/C]. Variation in the inferred [Si/C] with density 
depends on the UVB and rules out the quasar-only model unless [Si/C] increases sharply at low 
density. Comparisons with low-metallicity halo stars and nucleosynthetic yields suggest either that 
our fiducial UVB is too hard or that supermassive Population III stars might have to be included. 
The inferred [Si/C], if extrapolated to low density, corresponds to a contribution to the cosmic Si 
abundance of [Si/H]= —2.0, or flgi — 3.2 x 10~ 7 , a significant fraction of all Si production expected 
by z w 3. 

Subject headings: cosmology: miscellaneous — galaxies: formation — intergalactic medium — quasars: 
absorption lines 



1. INTRODUCTION 

Observational studies using high-resolution quasar ab- 
sorption spectra have firmly est ablished the presen ce of 
heavy elements such as carbon (jCowie et al.lll995T) . sili- 
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con l|Songaila fc Cowie Il99fil) . and oxygen l|Schave et alJ 
I2000a|) in the diffuse intergalactic medium (IGM). These 
metals constitute an important record of star formation 
and of the feedback of galactic matter into the IGM. 

This paper is the third in a series employing the 
statistics of pixel optical depths to study the enrich- 
me nt of the IGM. The b asic technique - pioneered 
by ICowie fc Songailal l)1998l see also Dav e et al. 1998 
and Songail a 199 8) and later employe d bv lEllison et alJ 
1)19991 l2000|) and iSchave et all l)2000a|) - was developed 
and extensively tested using cosmological hydrodynami- 
cal simulations by Aguirre, Schaye & Theuns (2001; here- 
after Paper I). The technique was then generalized and 
applied to 19 high-quality spectra in order to measure the 
full distribution of carbon as a function of redshift and 
gas density in S c haye et al. (2003; hereafter Paper II). 
See lAracil et"all pOOil and lPieri fc Haehnelti pOO^ for 
other recent studies applying the pixel method to C IV 
and O VI. 

Most studies of the enrichment of the IGM have fo- 
cused on C IV absorption because it is strong and lies 
redward of the Lya forest. Moreover, as shown in Paper 
II, ratios of C IV and H I optical depths can be converted 
into carbon abundances using an ionization correction 
that is neither very large, nor very sensitive to the tem- 
perature and density of the absorbing gas. While mea- 
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surements of the distribution of carbon provide impor- 
tant information on the mechanism by which the IGM 
was enriched, relative abundance information is crucial 
for identifying the types of sources responsible for the 
enrichment. 

Previous studies have established the presence of Si IV 
absorption in the IGM at z ~ 2 — 5 and have used simple 
ionization models to infer that their data is consistent 
with Si/C exceeding the solar ratio by a factor of a few 
ISongaila fc Cowie 11199ft ISongailalll998l ISongailall200ll 
Boksenberg, Sargent & Rauch 2003). The observed ra- 
tios of Si IV/C IV have also been used to study the shape 
and evolution of the ionizing U V background (U VB), but 
with conflicting results: while iSongailal l)1998f) sees an 
abrupt change in Si IV/C IV column density ratios at 
z ~ 3, Bokse nberg et all l|2003|) see no evidence for any 
evolution. 

This paper presents measurements of the relative abun- 
dances of silicon and carbon in the IGM, obtained by 
comparing the statistics of Si IV, C IV, and H I pixel op- 
tical depths in a sample of 19 high-quality quasar spectra 
to synthetic spectra obtained using a cosmological, hy- 
drodynamical simulation. Our primary goal is to mea- 
sure the overall Si/C abundance ratio in the gas for which 
Si IV absorption is detected, given a model for the extra- 
galactic ionizing UVB. In doing so, we also obtain some 
information on how much the distribution of silicon may 
differ from that of carbon (beyond an overall difference 
in the normalization), as well as constraints on the shape 
and the evolution of the UVB and on the thermal state 
of the absorbing gas. 

We have organized this paper as follows. In § 321 and [31 
we briefly describe our sample of QSO spectra and our 
methodology (both of which are discussed at length in 
Paper II). In S0]we first show results for our best QSO 
spectrum, Q 1422+230, as an illustration of the method. 
We then give measurements (using the full sample) of 
the Tsiiv/7"civ ratio, from which we infer the relative 
abundance of silicon to carbon ([Si/C]) for our fiducial 
UVB model (which is, as in Paper II, a re-normalized 
version of that given for galaxies and quasars by Haardt 
& Madau 2001, hereafter HM01). Next we give results 
for other UVB models, and discuss to what extent our 
data can constrain the UVB shape and evolution, and the 
variations in [Si/C]. We then give and interpret measure- 
ments of Tgav/THi and rsiin/rsav (which help constrain 
the thermal state of the absorbing gas). We discuss and 
interpret our results in SJ3 and conclude in [0 

2. OBSERVATIONS 

We analyze a sample of 19 high-quality (6.6 km s -1 ve- 
locity resolution, signal-to-noise ratio [S/N] > 40) ab- 
sorption spectra of 2.1 < z < 4.6 quasars. The sam- 
ple is identical to that of Paper II. It includes four- 
teen spectra t aken with the UV- Visu al Echelle Spectro- 
graph aiVES. IrrOdorico et all feOOOl on the Very Large 
Telescope (VLT) and five taken with the Hi gh Re solution 
Echelle Spectrograph (HIRES, IVogt et al.lll994[) on the 
Keck telescope. See Paper II (§ 2) for a description of 
the sample and data reduction. 

3. METHOD 

The basic method employed is similar to that described 
in Papers I and II, to which we refer the reader for de- 



tails and tests. Section 13.11 contains a brief overview of 
the method. Readers who want more detail would bene- 
fit from also reading sections 13.21 - 13.41 which summarize 
the method, making reference to the relevant sections of 
Papers I and II, and noting the changes to the method 
used in Paper II. 

3.1. Overview 

For each QSO spectrum, we first recover the optical 
depths due to H I Lya (A = 1216 A) absorption in all 
pixels in the Lya forest region. Then we recover the 
optical depth in pixels at the corresponding wavelengths 
of metal lines such as Si IV (AA1394, 1403), Si III (A1207), 
and C IV (AA1548, 1551) and correlate the (apparent) 
optical depth in one transition with that in another. We 
do this by binning the pixels in terms of the optical depth 
of the most easily detected transition, e.g., C IV, and 
plotting the median optical depth of the other transition, 
e.g., Si IV, against it. A correlation then indicates a 
detection of Si IV absorption; see the bottom curve in 
the left-hand panel of Fig. for an example. By doing 
the same for percentiles other than the 50th (i.e., the 
median) we obtain information on the full probability 
distribution of pixel optical depths (see the other curves 
in the panel). In this way a large quantity of information 
can be extracted from each observed spectrum. 

By comparing the observed optical depths with those 
obtained from synthetic spectra generated using a hy- 
drodynamical simulation, we make inferences about the 
distribution of silicon. For each observed spectrum, we 
generate a large set of spectra drawn from the simula- 
tion and process them to give them the same noise prop- 
erties, wavelength coverage, resolution, etc., as the ob- 
served spectrum. We do this for each of several UVB 
models; for each model the simulations include a carbon 
distribution as measured in Paper II for that assumed 
UVB, and some value of [Si/C], 10 with ionization bal- 
ances computed using the CLOUDY package 11 (version 
94; see Ferland et al. 1998 and Ferland 2000 for details). 
We analyze these simulated spectra in exactly the same 
manner as the observed ones, and compare the simulated 
and observed pixel optical depth statistics. 

The UVB models we employ are those in Paper II 
(§ 4.2). All are based on the models of Haardt & Madau 
(2001; see also 1996) 12 , but renormalized (by a redshift- 
dependent factor) so that the simulated spectra match 
our measurement (Paper II) of the evolution of the mean 
Lya absorption. Our fiducial model, "QG" , includes con- 
tributions from both galaxies (with a 10% escape frac- 
tion for ionizing photons) and quasars; "Q" includes only 
quasars; "QGS" is an artificially softened version of QG: 
its flux has been reduced by a factor of ten above 4 Ryd. 
Model "QGS3.2" is like model QG for z < 3.2 and like 
model QGS for z > 3.2 and was constructed to crudely 
model the possible evolution of the UVB if He II was 
suddenly reionized at z = 3.2. 

3.2. Recovery of pixel optical depths 

10 All abundances are given by number relative to hydrogen, and 
solar abun dance are taken to be (Si/H)© = —4.45 and (C/H)© = 
-3.45 I Anders fc OrevesselllM^ . 

11 See http : //www.nublado . org. 

12 The data and a description of the input parameters can be 
found at http://pitto.mib.infn.it/~haardt/refmodel.html. 
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After continuum fitting the spectra (Paper II, §5.1, 
step 1), we derive H I Lya optical depths thi for each 
pixel between the quasar's Lya and Ly/3 emission wave- 
lengths, although we exclude regions close to the quasar 
to avoid proximity effects (Paper II, §2). We use higher- 
order Lyman lines to estimate thi if Lya is saturated 
(i.e., F(X) < 3ct(A), where F and a are the flux and 
noise arrays, see Paper I, §4.1; Paper II, §5.1, step 2). 
For each H I pixel we then derive the corresponding Si IV 
(and C IV) optical depths rguv (and tciv), correcting for 
self-contamination and removing contamination by other 
lines (Paper I, §4.2). In addition, we recover Si III opti- 
cal depths rsiin, correcting for contamination by higher- 
order H I lines (Paper I, §4.2). We thus have four sets 
of corresponding pixel optical depths that may be corre- 
lated with each other. 

3.3. Generation of simulated spectra 

This study makes extensive comparisons between our 
set of observed spectra and those generated from a cos- 
mological simulation. The same simulation was used in 
Papers I (§3) and II (§4.1), to which the reader is referred 
for details. Briefly, the simulation uses a smoothed par- 
ticle hydrodynamics code to model the evolution of a 
periodic, cubic region of a (f2 m , ^a, ^6^ 2 , h, as, n, Y) — 
(0.3,0.7,0.019,0.65,0.9,1.0,0.24) universe of comoving 
size 12/i _1 Mpc, down to redshift z = 1.5 using 256 3 
particles for both the cold dark matter and the baryonic 
components. The UVB used in the simulation was cho- 
sen to match (and only affects) th e IGM temperature as 
measured bv lSchave et al.1 l)2000b|) . 

Synthetic spectra are generated as described in Papers 
I (§3) and II (§4.1), by computing the ionization balance 
for each particle using an assumed uniform UVB, then 
passing random sightlines through the snapshots of the 
simulation box and patching these sightlines together to 
form one long spectrum. Transitions due to C III, C IV, 
N V, Si III, Si IV, O VI, Fe II, and 31 Lyman lines of 
H I are included. The spectra also include noise, instru- 
mental broadening and pixelization chosen to match the 
observed spectra in detail. 

3.4. Comparison of simulations and observations 

Our si mulations (Paper I, §5.1) and tho se of oth- 
ers' (e.g.. lCroft et alJl997USchave et al.ll999|) show that 
there is a tight correlation between thi and the density 
and temperature of gas giving rise to the absorption. In 
Paper II we used this predicted correlation to compute 
an ionization correction (i.e., the ratios of C IV/C and 
H I/H) in order to convert ratios of tciv/thi into car- 
bon abundances as a function of density (Paper I, §6; 
Paper II, §5.1, step 5 and §5.2). This method works well 
for C IV because the ionization correction varies slowly 
with density (and therefore with thi) in the density range 
probed by the Lya forest (rin ~ 10~ 6 — 10~ 3 cm -3 ). Un- 
fortunately, this is not the case for Si IV: Fig. shows 
that while relatively insensitive to gas temperature, the 
ionization correction factor, [Si/H] — log(Tgiiv/Tm), in- 
creases dramatically with decreasing density, from < 10 2 
at an overdensity 6 = p/p ~ 10 (or njj ~ 10~ 4 cm~ 3 at 
z = 3) up to > 10 3 at 8 ~ 1. We find that this renders 
ionization corrections as used in Paper II unreliable for 
Si IV, because very small errors in the inferred H I optical 
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Fig. 1. — Ionization correction factor, [Si/H] — log(rgiiv/THi), 
as a function of temperature and hydrogen number density. Solid 
(dashed) contours are for the UVB model QG (Q) at z = 3. The 
hatched region indicates the temperature range containing 90% 
of the particles at the given density. For both backgrounds, the 
ionization correction increases rapidly for njj < 10 _4 cm — 3 . 



depth lead to large errors in the recovered silicon abun- 
dance. Rather than applying ionization corrections to 
the data, we therefore concentrate on directly comparing 
the predictions of our hydrodynamical simulation to the 
observed spectra. Using the measurements of the carbon 
distribution obtained in Paper II, we find the value of 
[Si/C] that best fits the data. 

For each observed quasar we generate 50 corresponding 
simulated spectra with the same noise properties, wave- 
length coverage, pixelization, instrumental broadening, 
and excluded regions. We apply to the simulated spec- 
tra the same automated continuum fitting routine as was 
used for the observations (Paper II, § 5.1, step 1; see also 
Paper I, § 4.1) to reduce any differences in the continuum 
fitting errors between simulated and observed spectra. 

Before generating the spectra, we first impose a car- 
bon distribution in the simulation that is consistent with 
the measurements presented in Paper II. In that paper 
we found that for a given overdensity and redshift, the 
probability distribution for the carbon abundance is well- 
described by a lognormal function, and we presented fits 
of the parameters of the distribution (i.e., the median 
and the width) as a function of overdensity and redshift 
for four different UVBs (Paper II, Eqn. 8 and Table 2). 
For the QG and Q ionizing backgrounds these fits are 
consistent with no evolution and we use the values for 
z = 3; for models QGS and QGS3.2 we use, for each 
quasar, the median redshift of the analyzed region. We 
divide each simulation snapshot into 10 3 cubes, and as- 
sign the gas particles in each section a carbon abundance 
of 

[C/H] = a + P(z - 3) + 7 (log<5 - 0.5) + s, 

where s, which is the same for all particles in the sub- 
volume, is drawn at random from a lognormal distribu- 
tion with mean and variance ct([C/H]) — 0.70, and 8 
is the overdensity of the particle. For QG, a = —3.47, 
7 = 0.65, and @ = Q; for Q, a = -2.91, 7 = 0.17, and 
(3 = (see Paper II, Table 2 for values used for other 
UVBs). Silicon abundances are assigned by assuming a 
constant value of [Si/C]. (The assumption that silicon 
tracks carbon perfectly is tested below.) 
For each observed and simulated quasar spectrum, op- 
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Fig. 2. — SilV optical depths in bins of tqiv (left panel) and thi (right panel) for Q1422+230. From bottom to top, the points are the 
median, 69th, 84th, 93rd and 97th percentiles. Solid lines represent predictions from simulations with ([C/H]) = — 3.8 + 0.65<5, <r = 0.70, 
and [Si/C]= log 5. For clarity, both observed points and simulated lines are plotted with vertical offsets of (from bottom to top) 0.0, 0.5, 
1.0, 1.5, and 2.0 dex. Note that for the right panel the spectra have been smoothed by 7.5 kms - 1 (see § 14.31 . 



tical depths for pixels corresponding to absorption by 
H I, C IV, Si IV and Si III are extracted as described 
in ii3.2l We may then plot various percentiles (such as 
the median) of the tsiiv distribution for each bin in tqtv- 
The same procedure can be applied to correlate Si IV 
with H I, or Si III with Si IV. For each C IV (or other) 
bin, errors on the observed points are calculated by boot- 
strap resampling the spectrum, i.e., we divide the Lya 
forest region of the spectrum into chunks of 5 A which 
are bootstrap resampled to form a new realization of the 
spectrum (see Paper II §5.1, step 3). We require that 
in each bin there be at least five pixels above the per- 
centile being computed, and that at least 25 pixels and 
five chunks contribute to the bin so that the errors are 
reliable; otherwise the point is discarded. For the sim- 
ulations we compute 50 synthetic spectra and each per- 
centile in each bin is set equal to the median of the 50 
realizations, with errors given by bootstrap resampling 
the 50 simulated spectra. We require each bin of each 
realization to have at least five pixels above the given 
percentile, and to have at least five pixels and one chunk 
contributing to the bin, and we discard medians com- 
puted with fewer than five acceptable realizations. 

For each percentile, the correlation disappears below 
some C IV optical depth r c at a value T m ; n that is deter- 
mined by noise, continuum fitting errors, and contam- 
ination by other lines. These may be corrected for by 
subtracting r m j n from the binned optical depths, thus 
converting points at thi < t c into upper limits (Paper II, 
section 5.1, step 4, and Fig. 4). For each realization, we 
compute T m ; n as the given percentile of the optical depth 
for pixels with thi < r c , where r c = 0.01 is chosen as the 
C IV or Si IV optical depth below which we never see 



a correlation. 1 The error on T m i n for each realization is 
computed by bootstrap-resampling the quasar spectrum. 
When the realizations are combined, the value of r m ; n is 
computed as the median among the realizations, and the 
error on this value is computed by bootstrap-resampling 
the realizations. 

The binned optical depth percentiles are then com- 
pared directly to the simulations, point by point. Be- 
cause of slight differences in contamination, accuracy of 
continuum fitting, and noise, T m - m for a given QSO may 
differ slightly from that of the simulations. Because we 
do not want to compare noise with noise, we either: A) 
subtract r m i n from both simulations and observations be- 
fore comparing, 14 or B) add a constant to the simulated 
optical depths, chosen to minimize the difference between 
the simulations and observations for thi < t c ; in this case 
we compare only points at thi > t c 

4. RESULTS 

Before presenting the results from our full sample, 
we illustrate the optical depth statistics by showing 
some results from our QSO with the strongest signal, 
Q1422+230. 

4.1. Results for Q 1422+230 

Fig-Elshows several percentiles of rsav; binned accord- 
ing to their corresponding optical depth in C IV (left) and 
H I (right) for Q1422+230. The bottom points are medi- 
ans, and the next four sets (with vertical offsets of 0.5,1.0, 

3 In Paper II we used functional fits to determine r c . For SilV 
the correlations are generally less strong than for CIV and we fix 
t c "by hand" . For bins in thi , we use t c = 1 . 

14 As in Paper II, we propagate the asymmetric errors for points 
in which T mln has been subtracted by computing a fine grid of 
errors (e.g., ±0.01cr, 0.02<r, . . .). 
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1.5, and 2.0 dex) represent percentiles 69.146, 84.134, 
93.319, and 97.725. These correspond to x=0, 0.5,1,1.5 
and 2cr values of the cumulative gaussian probability 
function f(x) = e _t / 2<T dt (so that, for exam- 

ple, a distribution of tsdv that is lognormal with width 
1 dex would give 1 dex higher metallicity in the 84.134th 
percentile than in the median). The solid lines are pre- 
dictions from the simulations with [Si/C] = log 5 = 0.70 
uniformly, median C metallicity ([C/H]) = —3.8 + 0.65(5 
and a width of cr([C/Hl) = 0.70 dex in the (lognormal) 
metallicity distribution. 15 

The left panel shows that the simulations with the 
QG UVB and [Si/C]= log 5 match the observations very 
well. This can be quantified by computing, for all points 
with logrciv > logr c = —2.0, the \ 2 difference between 
the simulated and observed data points. 16 We obtain 
X 2 /d.o.f . = 5.5/9, 4.0/9, 2.2/8, and 7.4/7 for the median, 

0. 5<7, lcr and 1.5<7 percentiles, respectively. 17 The good 
fit of the median indicates that a uniform [Si/C] ~ log 5 
model fits well (though slightly higher [Si/C] are favored; 
see below), and that the data do not require any strong 
trend of [Si/C] with density (for which rciv is a proxy). 
Comparison of the higher percentiles can help constrain 
scatter in tsuy/tciv, which depends on the scatter in 
both [Si/C] and the ionization correction, and it is clear 
that that simulations with uniform [Si/C] and a uniform 
UVB fit the observations well; this is quantified below 
using our full data sample. 

Comparison between simulated and observed rguy ver- 
sus thi (right) is somewhat less straightforward, as it de- 
pends both on [Si/C] and on the assumed distribution 
of carbon. Encouragingly, we find that for Q1422+230 
a model with uniform [Si/C] (and carbon distribution as 
derived in Paper II from the full sample) fits the observed 
medians fairly well. The predicted higher percentiles are, 
however, a bit low; this is partly because the Q1422+230 
C IV absorbers happen to exhibit slightly more metallic- 
ity scatter (« 0.81 dex; see Paper II, § 5.2) than the level 
derived from the full sample. 

In summary, the observed Q1422+230 Si IV optical 
depths, and their comparison to simulations show that: 

1. A uniform [Si/C]^ log 5 is favored, 2. There is no 
evidence for a trend of [Si/C] with density, 3. There is no 
evidence for scatter in rsnv / t~cw beyond that included in 
the simulation with uniform [Si/C] and a uniform UVB. 

4.2. rgiiv versus Tcrv f or the full sample 

To place more quantitative constraints on [Si/C], test 
for evolution, and make a more detailed comparison with 
several UVB models, we have combined the data points 
obtained from our entire sample. (These points are tab- 
ulated in Table 1.) Figure 0] shows log r siI v /tqtv versus 
logTcrVi m bins of z. To generate these points, we begin 
with Tgiiv values binned in tqw , as in Fig. [3 We then 
subtract from each the "flat level" r m i n for that QSO to 

15 These values are taken from the overall surface fits from Paper 
II; ([C/H]) is evaluated at z = 3 and the small redshift evolution 
is neglected; cr([C/H]) is evaluated at z = 3 and log S = 0.75, with 
both S- and ^-dependences neglected. 

16 We take into account both simulated and o bser ved errors, and 
adjust r m i n of the simulations as discussed in S 13.41 

17 The rather low reduced x 2 values may indicate some corre- 
lation between neighboring points, but they largely result from a 
slight overestimate of the errors near r c ; see Paper II, § 7. 
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Fig. 3. — Predicted optical depth ratios log Tsnv/7~CI v as a 
function of the temperature and the hydrogen number density, as- 
suming a solar ratio of Si/C. Solid (dashed) contours are for the 
UV-background model QG (Q) at z = 3. The hatched region indi- 
cates the temperature range containing 90% of the particles at the 
given density. 



adjust for noise, contamination etc. (see § EJ, then di- 
vide by the central value of the rciv bin. These points, 
gathered from all QSOs, are rebinned by x 2 -fitting, for 
each rciv range indicated in Fig.0] a constant level to all 
of the points in the specified redshift range. The errors 
represent 1- and 2-er confidence intervals (A% 2 = 2 and 
A X 2 = 4). 

The plotted lines indicate predictions from the simula- 
tions for our different UVBs; our fiducial model, QG, is 
shown in solid lines. For each background, we generate 
simulated Tsnv /^civ points in the same way as we did 
for the observations, but averaging over 50 simulated re- 
alizations as described in § 13.41 We then calculate a \ 2 
between all valid observed original (not rebinned) points 
and the corresponding simulated points. 18 Because we 
use 50 simulated realizations, the simulation errors are 
almost always negligible compared to the observed er- 
rors, but they are still taken into account by calculating 
the total x 2 using the formula: 



X = 



E 



Cobs 



where X = Tsav /t~cxv and a is the error in this quantity. 
We then add a constant offset to the simulated points 
(which corresponds to scaling [Si/C]) such that \ 2 is 
minimized. In each panel the lines connect the scaled, 
rebinned simulation points. 

The first clear result is that there is a strong trend of 
log(r S iiv/Tciv) with logr C iv, from < -2 at logr C iv ~ 
— 1.5 to > —1.0 at logTcrv —0.5. This previously un- 
noted correlation is expected if [Si/C] is constant and 
rciv correlates with density, since (as shown in Fig. |3J) 
Tsav/ r civ increases rapidly with gas density. As shown 
by the simulation lines, this trend is reproduced in all 
redshift bins by the simulations for all of our UVB 
models. Comparing the first three panels shows that 
while rsiiv/TCiv correlates strongly with tctv, there is 

18 Even using 50 realizations, it may occasionally happen that 
a simulated bins fails to have enough pixels for at least five re- 
alizations, and so is undefined; in this case the observed point is 
discarded as well. 
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Fig. 4. — Rebinned median log(rsirv/Tcrv) vs. logrciv in cuts of z for the combined QSO sample. The first three panels show bins 
centered at z = 2.0,3.0 and 4.0 with width Az = 1.; the bottom-right panel shows combined data from all redshifts. Lines represent 
corresponding rebinned simulation points (with errors suppressed, and with [Si/C] chosen to minimize the x 2 ) using different UVB models. 



Table 1. Recovered SiIV optical depths 



QSO 


log rciv 


lo 


5 TBUV 




0.5(7 


lo, 


? Tsnvi median 


logr SiIV ,+0.5cr 


lo 


S^SilV.+lo" 


lo 


;TSiiv, +i.5cr 


lo 


i-r sav ,+2a 


Q1422+230 


-9.000 










-3 


047 


0.088 


-2 


343 


0.035 


-1 


959 





031 


-1 


675 


0.035 


-1 


354 


0.047 


Q1422+230 


-3.900 


-3.371 





579 


-2 


530 


0.153 


-2 


115 


0.085 


-1 


898 





071 














Q1422+230 


-3.700 










-2 


827 


0.197 


-2 


183 


0.092 


-1 


817 





111 


-1 


622 


0.105 








Q1422+230 


-3.500 










-2 


847 


0.199 


-2 


308 


0.076 


-1 


978 





074 


-1 


765 


0.049 








Q1422+230 


-3.300 










-3 


214 


0.146 


-2 


492 


0.114 


-2 


054 





071 


-1 


786 


0.127 


-1 


310 


0.141 


Q1422+230 


-3.100 










-2 


879 


0.139 


-2 


270 


0.066 


-1 


918 





058 


-1 


605 


0.070 


-1 


326 


0.072 


Q1422+230 


-2.900 










-3 


007 


0.143 


-2 


361 


0.049 


-1 


929 





044 


-1 


583 


0.082 


-1 


293 


0.065 


Q1422+230 


-2.700 










-2 


936 


0.110 


-2 


303 


0.056 


-1 


954 





039 


-1 


729 


0.046 


-1 


482 


0.064 


Q1422+230 


-2.500 










-3 


076 


0.164 


-2 


306 


0.052 


-1 


936 





055 


-1 


661 


0.044 


-1 


390 


0.088 


Q1422+230 


-2.300 










-3 


127 


0.160 


-2 


371 


0.059 


-1 


985 





043 


-1 


648 


0.048 


-1 


402 


0.069 


Q1422+230 


-2.100 










-3 


194 


0.144 


-2 


409 


0.047 


-1 


966 





045 


-1 


691 


0.041 


-1 


344 


0.062 


Q1422+230 


-1.900 










-3 


050 


0.090 


-2 


315 


0.040 


-1 


974 





039 


-1 


697 


0.048 


-1 


450 


0.075 


Q1422+230 


-1.700 










-3 
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0.128 
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0.059 


-1 
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-1 
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0.050 


-1 


386 


0.063 
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-1.500 










-3 
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-2 
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0.077 


-1 


882 
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-1 
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0.094 


-1 


221 


0.121 


Q1422+230 


-1.300 










-2 


599 


0.173 


-2 


024 


0.125 


-1 


693 





120 


-1 


421 


0.104 


-1 


222 


0.086 


Q1422+230 


-1.100 


-3 


610 





606 


-2 


223 


0.321 


-1 


678 


0.213 


-1 


400 





165 


-1 


194 


0.143 








Q1422+230 


-0.900 


-3 


351 





772 


-1 


773 


0.439 


-1 


388 


0.193 


-1 


195 





121 


-1 


069 


0.118 








Q1422+230 


-0.700 


-1 


873 





313 


-1 


364 


0.203 


-1 


156 


0.102 


-1 


044 





067 


-0 


930 


0.095 








Q1422+230 


-0.500 


-1 


315 





497 


-1 


063 


0.270 


-0 


856 


0.226 


-0 


545 





225 














Q1422+230 


-0.300 


-1 


264 





813 


-0 


928 


0.450 


-0 


706 


0.325 























Q1422+230 -0.100 



NOTE. — Columns 1 and 2 contain the quasar name and the recovered CIV optical depth respectively. Columns 3 and 4 contain the 31st percentile of the recovered SiIV 
optical depth and the ltr error on this value. The remaining columns show the same quantities for higher percentiles (50th, 69th, 84th, 93th, and 98th). A CIV optical 
depth of log T^jy = —9 indicates that the corresponding TgilV are ^ ne T m in values. The complete version of this table, including all quasars, is in the electronic edition 
of the Journal. The printed edition contains only data for Q1422 + 230. 
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Fig. 5. — Rebinned median log(rgiiv/Tciv) vs - z m cuts of rciv for the combined QSO sample. The first three panels show bins centered 
at logTciv = —1.25,-0.75 and -0.25 with width 0.5 dex; the bottom-right panel shows data for all tqiv combined. Lines represent 
corresponding rebinned simulation points (with errors suppressed, and with [Si/C] chosen to minimize the x 2 ) using different UVB models. 
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Table 2. Best fit [Si/C] and x 2 /d.o.f. 



UVB model best fit [Si/C] x 2 /d-°-f- 



QG 


n 77 +0.05 
u -' '-0.05 


65.7/115 


Q 


1 -^°-0.06 


65.6/115 


QGS 


26+ ' 06 
u ' ZD -0.07 


73.8/115 


QGS3.2 


046 +o.io 


85.2/115 



little dependence on redshift. This can be seen more 
clearly in the first three panels of Fig. [5] which show 
log(Tsiiv/T"civ) versus z in bins of tciv- There is no ev- 
idence, in either the simulated or the observed points, 
for evolution in tsuv/tciv, except perhaps for a slight 
increase in rsiiv/ T civ with increasing z at the highest 
densities. 

The observed trends in log(rsiiv/Tciv) are reproduced 
well by the simulations. Because Tsiiv/ T 'civ scales with 
Si/C, the offset in tsuy/tcw obtained by minimizing 
the x 2 (Eq. against the observations can be used 
to compute the best fit [Si/C]. For our fiducial UVB 
model QG, the simulated spectra were generated with 
[Si/C]=0.70, we find an offset of +0.07 dex (implying 
a best-fit [Si/C]=0.77) with x 2 /d.o.f.=65. 7/115. As we 
found in Paper II and for Q1422+230 above, the reduced 
X 2 is somewhat low; in Paper II we showed that his was 
largely due to a slight overestimate of the errors at low- 
tciv- 

The fitted [Si/C] values and corresponding x 2 /d.o.f., 
are listed in Table|21 with errors computed by bootstrap- 
resampling the quasars used in the x 2 minimization. For 
our fiducial model, QG, the best fit [Si/C]= 0.77±g gg. 
The quasar-only background Q (which is probably too 
hard; see Paper II) gives a higher values of [Si/C] = 
1.48±g;g|, and the softer QGS and QGS3.2 backgrounds 
give lower values than QG by w 0.3 — 0.5 dex. Note that 
the QGS background is unrealistically soft at z ^ 3 (see 
Paper II), but the QGS3.2 background may be plausible. 

Because the Tsiiv/T'civ ratio depends on both [Si/C] 
and the shape of the UVB, its measurement can be 
used to study the evolution of th e UVB under th e as- 
sumption that [Si/C] is constant. ISongailal l)1998L see 
also Songaila & Cowie 1996) has measured the median 
Si IV/C IV ratio versus z for C IV systems of col- 
umn density 5 x 10 12 c m - 2 < V(CIV) < 10 14 c m - 2 
and found evidence for strong evolution, as well as a 
sharp break in tsuv/tctv &t z — 3.0. This was inter- 
preted as evidence for a sudden soft ening of the UVB 
at z > 3. However, in a recent study IrToksenberg et alJ 
(2003) , in agreement with their earl ier work and that 
of lKim. Cristiani. &: D'Odoricol l)2002|) find no evolution 
in rsiiv/7"civ in their sample of 10 12 cm~ 2 < AT (CIV) < 
3 x 10 14 cm~ 2 absorbers. As discussed above (see Fig.[5J), 
we see no evidence for evolution in the median rsnv /tciv 
stronger than that predicted by the simulations, in any 
of our tciv bins, or when all tciv values are combined 
(Fig. El bottom-right). 

Because rsiiv/ r civ varies by ~ 1.5 dex in correlation 
with tciv j hs evolution is best assessed by using only a 
small window in tciv; otherwise evolution in the weight 
provided by each tciv ~~ whether due to selection effects 



or evolution in the distribution of tciv ~~ may lead to ap- 
parent evolution in Tsiiv/T'civ- I n our analysis, in which 
we use small cuts in tciv > 19 we see no evidence for evo- 
lution in the UVB beyond that in the smoothly changing 
QG and Q models; indeed, the QGS3.2 model with an 
artificial change in softness at z = 3.2 is disfavored by 
our data, having a significantly higher x 2 than either the 
QG or Q models. 

The simulations we employ assume a constant and uni- 
form [Si/C]; but because the formation mechanism for 
Si and C may be different this need not be the case. 
We can, however, observationally constrain the scatter 
in [Si/C] by repeating our determination of it using dif- 
ferent percentiles in tsuv/tcyv- If the probability dis- 
tribution of [Si/C] for fixed tciv is> like that of [C/H] 
and [Si/H] for fixed thi, lognormal, then we can di- 
rectly constrain the width <r([Si/C]) of the distribution 
by comparing [Si/C] derived for different percentiles. For 
the 69th, 84th, 93rd and 97th percentiles, we obtain 
[Si/G]=0.73±°;gg, 0.64t°$, 0J58l°;& and 0.781°;!°, re- 
spectively. This translates into a rough 2a upper limit 
of <r([Si/C]) < 0.04 dex. 20 As a rough test of this upper 
limit, we have generated simulated spectra with a median 
[Si/C] = 0.77 and the usual C distribution, but with a 
lognormal scatter in [Si/C] of width <r( [Si/C]). If we add 
the x 2 for the 69th, 84th and 93rd percentiles for quasars 
at z > 3 (there is little useful information on the upper 
percentiles from z < 3), we find x 2 /d.o.f. = 198.0/185 
for no scatter in [Si/C], and x 2 = 198.2, 203.4, and 223.0 
for cr([Si/C]) = 0.1,0.2,and 0.4 dex, respectively. Be- 
cause the percentiles are correlated these cannot be cor- 
rectly translated directly into confidence limits; how- 
ever they suggest that the data are compatible with 
cr([Si/C]) = 0.1, but probably not with <r([Si/C]) = 0.2 
and almost certainly not with <r([Si/C]) = 0.4. Thus, 
contrary to the large scatter in [C/H], <r([C/H]) ~ 0.7dex 
found in Paper II, there appears to be very little scatter 
in [Si/C]. 

We may also subdivide our sample by redshift and den- 
sity to test the dependence of [Si/C] on these. First, 
computing [Si/C] using only spectra that have a median 
absorption redshift med(2;) > 3.0 (see Paper II, Table 
1) yields [Si/C]=0.76±8:o7, versus [Si/C]=0.79t°;J° us- 
ing the spectra with med(z) < 3.0. The [Si/C] values in- 
ferred from the redshift subsamples are also (marginally) 
consistent for the Q and QGS UVBs, but inconsistent 
for QGS3.2; the latter would imply a jump from [Si/C] = 
0.77t°,: 9 at z > 3 to [Si/C]= 0.32t°;°j at z < 3. This 
is a second way of seeing that our data disfavors a sud- 
den change in UVB hardness near z — 3, assuming that 
[Si/C] is constant. 

To test for variation in [Si/C] with overdensity S we 
have recomputed [Si/C] using only pixels with thi corre- 
sponding to 5 < 20 or 6 > 20 (using the thi-6 conversion 
of Paper II, Fig. 2). We obtain [Si/C]=0.87lg;Jg for 
the low-density sample, versus [Si/C]=0.73t°;os at high 

19 In the bottom right panel of Fig. |3 we combine all tciv 
values. Although the observations and simulations should have 
similar weightings by C IV, some differences may remain and the 
comparison is less reliable than if cuts in tciv are made. 

20 This is only a rough error estimate as it assumes that the 
points are fully independent, which does not hold because the 
points correspond to a ranking. 
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= 0.25 

= 0.48 

1.0 1.5 2.0 
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Fig. 6. — Median log(rsiiv/THl) vs. logTfn in bins of z for the combined QSO sample. The first three panels show bins centered at 
z = 2, 3, and 4 with width Az = 1; the bottom-right panel shows combined data for all redshifts. The lines represent corresponding 
simulation points (with errors suppressed) using different UVB models. 



density. This difference could imply either that [Si/C] is 
higher at low densities, or that the UVB is softer than 
we have assumed (see Fig. |3J; but the effect is only sig- 
nificant at the 1.3a level. 

4.3. rsiiv versus thi for the full sample 

While the rsiiv/^civ ratios give the most direct con- 
straints on [Si/C], it is also useful to examine tsuv/thi: 
comparing the simulated to the observed tsuv/thi ratios 
gives, in principle, a second estimate of [Si/C]. This infer- 
ence, however, is less reliable than that from ts;iv/ t civ 
for two reasons. First, it includes the uncertainty in 
[C/H], which is largest at the relatively high densities 
at which (because of the small Si IV/Si fraction at lower 
densities) Si IV is best-detected. Second, Si absorption 
lines have a much smaller thermal width than hydrogen 
lines, and are only observable in relatively high-density 
gas (where Hubble broadening is small). This leads to 
significant differences between the Si IV-weighted and 
H I-weighted densities (see Paper II, § 4.3, Paper I, § 5.3), 
i.e., the Si IV and H I absorption do note arise from ex- 
actly the same gas. The effect of this is to skew (and 
weaken) the correlation of rgnv with thi in a way that 
depends on the H I colu mn-density distributio n. Unfor- 
tunately, as discussed bv lTheuns et alJ C2002), the sim- 
ulation does not exactly reproduce the observed iV(HI) 
distribution at the high-column density end (which is 
most important for the present study), so the "skewing" 
may be somewhat different in observed and simulated 
correlations, and render conclusions about [Si/H] unreli- 



able. 

The problems arising from differential line broaden- 
ing can be partially remedied by smoothing the spectra 
so that the minimal line widths of all species become 
similar; experimentation shows that smoothing the spec- 
tra by convolving them with a Gaussian with FWHM of 
~ 5 — 10 km s" 1 significantly increases the strength of the 
signal, and a smoothing of 7.5 kms -1 has been adopted 
for the calculations shown in Figs. [5] and 03 Fig. shows 
log rsiiv / t hi versus log thi in bins of z for our combined 
sample. Lines again connect the corresponding simula- 
tion points (with an overall scaling to best match the ob- 
servations) that reproduce the observed trends in z and 
thi- The scalings correspond to best- fit [Si/C] values of 
0.97 ± 0.08, 1.65 ± 0.08, 0.25 ± 0.07 and 0.48 ± 0.09 for 
QG, Q, QGS, and QGS3.2 respectively, slightly higher 
than the values found above using tsuv/tcw- However, 
as explained above, these values depend on the degree of 
smoothing: for example, with no smoothing, we recover 
values sa 0.1 — 0.3 dex higher yet. Since it is unclear which 
smoothing level gives the correct results, the inferences 
of [Si/C] using Tsnv/THi should not be relied on. 

The same difference in thermal width affects inferences 
using Si IV/C IV, but at a lower level; whether or not 
we smooth by 7.5 kms" 1 changes the inferred [Si/C] by 
~ 0.1 dex, which is a reasonable estimate of the induced 
systematic error. 

4.4. Tgiiii versus Tsnv for the full sample 
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Fig. 7. — Predicted optical depth ratios log Tsiiu/rgiiv as a func- 
tion of the temperature and the hydrogen number density. Solid 
(dashed) contours are for the UV-background model QG (Q) at 
z = 3. The hatched region indicates the temperature range con- 
taining 90% of the particles at the given density. 



A final correlation that we have examined is that 
of Tsini/Tsiiv with rsiiv- Fig- 13 shows the predicted 
T'Siiii/T'siiv versus the temperature T and density bh of 
the absorbing gas. For high T, Tsnn/rsiiv becomes in- 
dependent of rin, and declines rapidly with T. Thus, the 
presence of Si III can be used to constrain the tempera- 
ture of the gas providing Si IV absorption. In Paper II 
a similar test was carried out using the tciii/tciv ratio, 
yielding the constraint T < 10 5 K. 

Fig. [S] shows Tsaii/rsiiv versus rsiiv for our full sam- 
ple, for three cuts in z. At logrsuv ^ —1.2 and 
2.5 < z < 3.0 (where the data are particularly good) 
—0.5 < Tgini/Tsiiv ^ 0.5, corresponding to a direct up- 
per limit of T < 10 4 - 9 - 10 5 1 K for the bulk of the gas 
giving rise to this Si IV absorption (the lower value per- 
taining to the higher end of the rsnv range.) The sim- 
ulations give predictions that depend only very weakly 
on the UVB model and are in good agreement with the 
data, particularly at 2.5 < z < 3.0 (however note that 
at z > 3.0 there is somewhat more Si III absorption pre- 
dicted than observed for high rsuv)- Concentrating on 
the QG model, we find x 2 /d.o.f. = 80.5/114 comparing 
all simulated and observed points with logrsiiv > —2.0. 
Fitting an offset to the median logTsnii/rsiiv ratios (as 
done above for rciv/ T Siiv) we find — 0.32+Q'![g dex, in- 
dicating that simulations predict somewhat too much 
Si III absorption overall. We have repeated this exer- 
cise for the 16th, 31st, 69th and 84th percentiles in order 
to fit the center and width tr([SiIII/SiIV]) of a lognor- 
mal distribution governing scatter in Tsim/rsiiv beyond 
that present in the simulations. The fit obtained indi- 
cates that overall the observed Tsiin/Tsiiv ratio is lower 
by 0.15 ± 0.06 dex than that for the simulations, with 
cr([SiIII/SiIV]) = 0.07 ± 0.07 dex. These results imply 
that most of the gas is cool enough to be consistent with 
the photoionization equilibrium assumed in the simula- 
tions, and that the scatter in rsan/rsiiv is similar to that 
in the simulations. But they are also consistent with a 
(small) contribution by hotter gas to the observed Si IV 
optical depths. 

5. ANALYSIS AND DISCUSSION OF RESULTS 
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Fig. 8. — Ratios of median Tsini/Tgiv versus rsuv in bins of z 
for the full sample. The lines represent corresponding simulation 
points (with errors suppressed) using different UVB models. 



We have shown that the pixel optical depth correla- 
tions derived from our observed QSO spectra are consis- 
tent, in detail, with spectra drawn from a hydrodynam- 
ical simulation with: 1) an assumed carbon metallicity 
distribution as derived in Paper II from measurements 
of C IV abs orption, 2) a uniform (re scaled) UVB model 
taken from Haardt & Madau (2001), and 3) a constant 
and uniform [Si/C] value. 

We have used this success to draw inferences regarding 
[Si/C] for each of several models of the UVB. Before ana- 
lyzing these inferences it is worth discussing some effects 
that bear on their robustness. 

5.1. Uncertainties 

UVB spectral shape: Our inference of [Si/C] obviously 
depends on the assumed shape and evolution of the UVB: 
the models used here produce a range 0.26 < [Si/C] < 
1.48. These models, however, span a range of possi- 
bilities that is larger than that allowed by independent 
observations. The harder model Q, for example, was 
found in Paper II to produce a mean carbon abundance 
that increases with decreasing density, which is probably 
unphysical. Likewise the softer model QGS was found 
to imply increasing [C/H] with z, also unphysical. The 
model QGS3.2, with a sudden transition at z = 3.2, pro- 
duced no problems in Paper II, but here we find that 
the resulting predicted jump in tsiiv/tciv is disfavored 
by our observations. Our fiducial model QG appears 
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compatible (and produces physical results) for all of our 
metallicity measurements thus far; it therefore seems un- 
likel y that this UVB is very f ar in error. Note, however, 
that iBoksenberg et al"l l|2003j) found that they were un- 
able to fit Si II / Si IV vs. C II/C IV column density ratios 
in their data (probing higher densities than our analysis) 
using Haardt & Madau (2001) models. 

UVB inhomogeneity: We have assumed that the UVB 
is perfectly uniform, which is probably unrealistic. How- 
ever, we can limit the non-uniformity of the hardness 
of the UVB because such non-uniformity would lead to 
scatter in rsiiv/ r civ even for uniform [Si/C], whereas 
our measurement of the higher percentiles in tsuv /''"civ 
proved to be compatible with the predictions of the sim- 
ulations for a uniform [Si/C]. 

Collisionally ionized gas: The simulation we have em- 
ployed does not contain feedback from star formation and 
therefore contains only the hot (T ^> 10 5 K) gas result- 
ing from heating by accretion shocks. If the true universe 
has hot gas not present in our simulations, our inferences 
could be biased. We have put a constraint T < 10 4 9 K on 
the temperature of the bulk of the gas giving rise to the 
Si IV absorption, using the rsnn/Tsav ratio. The simu- 
lation, however, does predict slightly higher rsim/Tsiiv 
overall, by ps 0.3 dex using the median and by ks 0.15 dex 
using all percentiles. Although this might be attributed 
to some systematic difference between simulations and 
observations, it may indicate that some fraction (< 30%, 
since that is the lowest percentile we probe) of the Si IV 
absorbing gas is hotter (T > 10 5 K) than in the sim- 
ulations. Note also that, as in Paper II, we can put 
no constraint on metals in gas that is either too hot 
(T > 10 5 K) or too cold (T < 10 4 K) to cause detectable 
Si IV or C IV absorption. 

Uncertainties in recombination rates: iSavinl l)2000fl 
has pointed out that the Si IV and C IV dielectric 
recombination rates used in CLOUDY may be uncer- 
tain by up to ~ 0.3 dex, leading to uncertainties in 
inferred [Si/Cj. Conside r ing fi rst C IV, the laboratory 
work of Schip pers et alJ l)2001lL with quoted uncertain- 
ties of ~ 15%, finds rates ~ 15 — 30% higher than those 
used i n CLOUDY version 94. The calculation of Savin 
( 2000) suggest that if the CLOUDY rate were 30% too 
low, we would overestimate [Si/C] by a factor that de- 
pends on the density of the absorbing gas and varies from 
~ 0.0 — 0.07 dex for the range given below in § 15.21 For 
Si IV the situation is worse because there are no pub- 
lished laboratory ex periments, s o the rate could be un- 
certain by ±0.3 dex ijSavinl 120001. If the CLOUDY rate 
were too high (resp. too low) by 0.3 dex, we would over- 
estimate (resp. underestimate) [Si/C] by 0.2 dex at the 
highest gas densities we probe, by a negligible amount at 
the lowest densities, and probably by ~ 0.1 dex overall. 

Note also that differences in the Si and C line widths, 
as discussed in § 14.31 could lead to an uncertainty of ~ 
0.1 dex in inferred [Si/C] using our method. 

5.2. Corresponding gas and column densities 

We have analyzed Si IV, C IV, and H I pixel optical 
depths to draw inferences regarding the cosmic silicon 
abundance. For these to be compared with those from 
theoretical or observational studies, it is useful to esti- 
mate the gas densities and the H I and C IV line column 
densities our measurements pertain to. 



In measuring rsnv/Tm vs. thi, (Fig. EJ, we obtain 2a 
detections for 0.75 < logrni < 2.5 for 2.5 < z < 3.5 
(where our data are best). Using the theoretical thi- 
density relation given in Paper II, this range corresponds 
at z = 3 to 0.5 < log 8 < 2.0. (Note, however, that the 
upper limit results from the inability to recover high- 
thi pixels and does not apply to results binned in rciv-) 
These densities can be converted into approximate H I 
column densities using the theoretical relation iVm ~ 
10 15 cm- 2 (<S/10) 3/2 [(l + z)/4] 9 / 2 from lSchavelll200l . to 
obtain a column density range 

2 x 10 14 cm" 2 < N m < 3 x 10 16 cm" 2 . (2) 

Another way to estimate the H I column densities cor- 
responding to the H I optical depth range over which we 
detect Si IV is to use the relation between column density 
TV and optical depth r c at the center of a line of width b, 

Tc = °' 76 OfflAo.Hi) ( 10" cm- 2 ) (l0kms-i) ' 

(3) 

where / and Ao are the oscillator strength and the rest 
wavelength respectively. Using b = 25 kms -1 yields 
2 x 10 14 cm" 2 < Nm J$ 1 x 10 16 cm -2 , in close agree- 
ment with our earlier estimate. We can estimate the 
C IV column density range in a similar way. From 
Fig. 0] we see that Si IV is detected over the range 
— 1.5 < logrciv 0-0, which corresponds to 

7 x 10 n cnV 2 < Vciv < 2 x 10 13 cm" 2 , (4) 

for a typical C IV line width of 10 kms -1 . 

Note also that our analysis pertains to signif- 
icantly lower physical densities than results ob- 
tained using the C II/C IV ratio as a density 

tracer llSongailal Il998t IBoksenberg et all [2003); using 
Fig. 23 of IBoksenberg et alJ l)2003j) . our upper density 
limit corresponds roughly to V(CII)/V(CIV) < -1.5. 2 ' 1 

5.3. Comparison to Nucleo synthetic yield studies 

The small scatter in [Si/C] implies that silicon and car- 
bon share a common enrichment mechanism, and both 
the high [Si/C] and early enrichment epoch we observe 
suggest enrichment by massive stars. It is therefore inter- 
esting to compare our [Si/C] to that predicted in theo- 
retical supernova yields for low-metallicity, zero metal- 
licity, and (more speculatively) supermass ive progeni- 
tors. The yields ofWo oslev fe Weaver! l)1995j) . integrated 
over a Salpeter initial mass function (IMF) with slope 
between -2.5 and -0.5, give production factors (yields 
relative to solar metallicity) of 0.4 < [Si/C] < 0.6 for 
10~ 4 < Z/Zq < 0.1 progenitor met allicities. At Z = 
signifi cantly smaller [Si/C] results: IWooslev fc Weaver! 
lll995h and IHeger fc WoosTevI ( 2002) give production 
factors [Si/C]~ 0.15 for massive 12-40 M Q stars, and 
[Si/C]- 0.0 - 0.2 for 15 - 30 Mm met al-free stars are 
also quoted bv lChiefn fc Lim onfiil l|2flQ2| ) from bo t h thei r 
own calculations and those of lUmeda fc Nomotol (2002) . 
Zero metallicity also, however, allows (at least theoreti- 
cally) for supermassive s tars with quite different yields. 
IHeger fc WoosTevI (120751 provide yields of [Si/C]~ 1.2 

21 I f the relatively high-density gas probed by Bokscnbcr g et all 
12003) (and Songaila 1998) lies relatively close to galaxies, it may 
be strongly affected by those galaxies' soft UV radiation. 
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for supermassive (140-260 M Q ) progenitors that produce 
pair-instability supernovae. Including as well the mas- 
sive (12-40 M ) stars with a Salpeter IMF gives [Si/C]~ 
0.85. 

We may also compare our [Si/C] values to those in 
metal-poor stars, but here there are a range of results, 
and substantial uncertainties - particularly with regard 
to C - concerning the impo rtance of mixing and accretion 
from a companion (e.g., iCarretta. Gratton. fc Snedenl 
I2000j) . The lowest-metallicity ([Fe/H] < -2.5) 

halo stars give typically [Si/C]~ —0.4 1-0.5 (e.g., 

Chieffi fc Limongii 20021 iNorris. Ryan, fc Beerl 120011: 
Denagne et alJ2002HMc William et alJl995lh For -2.5 < 
[Fe/H] < —1.5, somewhat higher values are appar- 
ently seen: measurements of —0.5 < [C/Fe] < 0.0 and 
0.2 < [Si/Fe] < 0.7 (see the compila tion by Norris et 

ah 2001) suggest [Si/C] - 0.7, and IMc William eTaD 

lITQQl finds values of [Si/C] ~ 0.2 - 1.2. How- 
ever. ICarretta. Gratton. fc Snedenl l)200f)fl find higher 
[C/Fe] of —0.2 < [C/Fe] < 0.3 when using only unmixed 
stars, implying somewhat smaller [Si/C] ~ 0.4. 

If we assume that our fiducial UVB model (and as- 
sumed recombination rates; see > I5.1JI are correct, then 
the inferred [Si/C]= 0.77 ± 0.05 has interesting impli- 
cations. It is substantially higher than the [Si/C] pre- 
dicted from < 40 M Population III stars with a stan- 
dard IMF, and it is somewhat (~ 0.2 dex) higher than 
abundance ratios in low-metallicity stars, and predicted 
yields for massive Population II stars. If the latter dis- 
crepancy is taken seriously, it could be remedied by in- 
cluding some s upermassive (pair-instab ility) supernovae 
as modeled bv lHeeer fc Wooslevl (|2002i). 22 

If, instead, we assume that [Si/C] should be ~ 0.5 on 
the basis of standard nucleosynthesis yields, then we infer 
that the UVB must be softer (at bith low- and high-z) 
than our fiducial model QG. However, it should be noted 
that a UVB as soft as model QGS, for which we find 
[Si/C]= 0.26, was demonstrated to produce unphysical 
results in Paper II. 

5.4. Implications for cosmic abundances 

In §01 we measured [Si/C] values for several UVB mod- 
els by comparing observations of tsuv/tcw to simula- 
tions generated using a given UVB model and a distri- 
bution of [C/H] as determined in Paper II (§ 7) for that 
UVB. Combining these fits (See Table |2J with the re- 
sults of Paper II enables us to assess the overall cosmic 
abundance of silicon. For these calculations we assume 
that our derived [Si/C] is uniform and constant, but we 
include the uncertainty in its determination. 

In Paper II we combined the median [C/H] (5, z) with 
the width er([C/H])(<5, z) of the lognormal probability dis- 
tribution of [C/H] for —0.5 < log S < 2.0 to determine 
the mean C abundance versus S. This was then inte- 
grated over the mass- weighted probability distribution S 
(obtained from our hydrodynamical simulation) to com- 
pute the contribution by gas in this density range to the 
overall mean cosmic [C/H]. Assuming that [Si/C] is con- 
stant over this density range 23 we obtain, for our fiducial 

22 That supermassive stars might be necessary to explain our 
observed [Si/C] value is particularly interesting in light of the ap- 
preciable fraction of all cosmic metal pro ducti on at z > 3 that may 
be inferred to lie in the Lyct forest (see § 15.41 . 

23 Note that this is an extrapolation beyond the range over which 



UVB model QG, [Si/H]= -2.03 ±0.14, corresponding to 

n Si ,iGM - 3.4 x icr^io™ 2 - (J^j , (5) 

with no evidence for evolution. Extrapolating our [C/H] 
and [Si/C] results to the full density range of the simu- 
lation would yield values « 0.2 dex higher. For a harder 
UVB, the inferred cosmic Si abundance would be sig- 
nificantly higher, as both the inferred [C/H] and [Si/C] 
increase with the UVB hardness. For our quasar-only 
model Q, we obtain a contribution o f [Si/Hjsa —0.8, or 
Q Si w 5 x 10~ 6 . For comparison, Songaila (200J) 
found, using a direct sum of Si IV lines, that f^siiv — 
(0.5 — 2) x 10~ 8 . This indicates that although Songaila's 
results provide a useful lower limit, an ionization cor- 
rection of more than 10 (or even much larger for harder 
backgrounds) must be made to translate f^siiv into an 
estimate of the true cosmic silicon abundance. The large 
ionization correction implies that the evolution of fisiiv 
does not by itself provide interesting constraints on the 
evolution of fisi- 

It is interesting to compare our derived metallicities 
to the metal production expected of high-z galaxies and 
contained in other known types of objects. Integration 
of the observed cosmic star formation rate indicates that 
fi*(z = 2.5) s» 0.001, a quantity sufficient to supply a 
contr ibution of 1/30 solar metallicity (see, e.g.. iPettini I 
2003). Thus if Si is used as a tracer of metals, then 
the metals observed in the Lya forest would represent 
w 30 — 50% (for the QG UVB, depending on whether we 
extrapolate to high-densities) of all metals produced by 
z = 2.5. If the Q UVB were to hold, the IGM would hold 
~ 5 — 8 times the expected metal production, which is 
probably yet another indication that the Q UVB is unre- 
alistically hard. Another, similar approach is to compare 
our derived f^si to the total amount of Si locked in stars: 

where Zgi is the mean silicon metallicit y in those stars , 
and f2» w 0.001 is the estimate given bv iPettinil (2003) 
for the stellar density at z = 2.5. Thus, even if all stars 
at z « 2.5 were to carry solar abundance of silicon, they 
would contain - 60 - 70% (for QG) or - 8 - 12% (for 
Q) of the cosmic silicon; the rest would be stored in the 
IGM. 

It is also worth noting that, if silicon is used as the 
metallicity tracer and our fiducial UVB is adopted, the 
Lya forest contains about twice the cosmic metal mass 
pr ovided by dam ped Lyman absorbers, as estimated 
bv IPettini 1 12003]) . 

6. CONCLUSIONS 

We have studied the relative abundance of silicon in 
the IGM by analyzing Si IV, C IV, and H I pixel opti- 
cal depths derived from a set of high-quality VLT and 
Keck spectra of 19 QSOs at 1.5 < z < 4.5, and com- 
paring them to realistic, synthetic spectra drawn from 
a hydrodynamical simulation to which metals have been 

we have measured [Si/C]. 

24 As in Paper II, we have corrected Songaila's value, which as- 
sumes an Q m = 1.0 cosmology to the Q m = 0.3 cosmology adopted 
here. 
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added. Our fiducial model employs the ionizing back- 
ground model ("QG") taken from Haardt & Madau 
(2001) for quasars and galaxies (rescaled to reproduce 
the observed mean Lya absorption) , and assumes a car- 
bon abundance as derived in Paper II: at a given density 
5 and redshift z, [C/H] has a lognormal probability dis- 
tribution centered on —3.47 + 0.65(log£ — 0.5) and of 
width 0.70 dex. The main conclusions of this analysis 
are as follows: 

• For our fiducial model, the median optical depth 
ratios Tsiiv/Tciv for —1.5 < logrcrv are repro- 
duced well as a function of z and tqw by the sim- 
ulations, if we take [Si/C]=0.77± 0.05, uniformly 
and at all times. These measurement pertain to gas 
(over)densities of S > 3, or column densities -/Vhi > 
2 x 10 14 cm^ 2 and AT C rv ~7x 10 11 - 2 x 10 13 cm" 2 
at z = 3. 

• We find a strong correlation between tsuv/tcw 
and tciv exhibited by the data that are reproduced 
by the simulations. This indicates that evolution 
in Tsnv/Tciv is best assessed using small cuts in 
tciv- Our results agree very well with simulations 
using a smoothly evolving UVB and show no sig- 
nificant evolution in tsuv/tcw (except perhaps a 
slight rise in tsuv/tcw with z for logrciv ^ —0.5). 
T he lack of evolution is co nsistent with the results 
of iBoksenberg et all ll2003ft and Kim et al. (2002), 
but not those of lSongailal l)1998|) . 

• The [Si/C] value inferred from rsnv/Tciv depends 
on the assumed UVB. A harder UVB model Q (us- 
ing only the contribution of quasars to the UVB) 
gives a much higher ratio, [Si/C]=1.48loQg. A 
model "QGS3.2" in which the UVB is decreased 
in intensity by 1 dex blueward of 4 Ryd at z > 3.2 
(crudely simulating incomplete He II reionization) 
gives a lower ratio 0.461°,;^ (or 0.321°;°* using the 
z > 3 data alone), but gives too much evolution in 

TSilv/lCIV- 

• Subdividing the sample by z and gas density 5, 
we see no evidence for evolution in [Si/C] for any 
of our UVBs (except QGS3.2, which gives signif- 
icantly different results at low- and high-z). The 
low-<5 and high-<5 [Si/C] inferences are discrepant 
by 1.4(7 for our fiducial UVB and by 3.8tr for the 
quasar-only UVB. This indicates either that [Si/C] 
increases with decreasing density, or that our fidu- 
cial UVB may be slightly too hard (and that the 
QSO-only UVB is far too hard). 

• While the dominant uncertainty in [Si/C] comes 
from the shape of the UVB, an additional uncer- 
tainty of up to ~ ±0.1 dex in our inferences (and 
those of other studies) could result from uncertain- 
ties in the assumed Si IV dielectric recombination 
rates. There may also be a systematic uncertainty 
of up to ~ 0.1 dex resulting from the different ther- 
mal widths of C and Si. 

• Evaluating [Si/C] for optical depth percentiles 
higher than the median gives no evidence for addi- 
tional scatter, either in Si/C or in the hardness of 



the UV background (and hence the ionization cor- 
rection). The width of a lognormal distribution 25 
of [Si/C] is constrained to be much smaller than 
that of [C/H]. 

• Analysis of tsuy/tui also gives information on 
[Si/H] and [Si/C] but is subject to large systematic 
effects because of the significantly different thermal 
width of Si and H. If the data is smoothed to min- 
imize this difference, we can roughly reproduce the 
inferences based on rsnv/Tciv- 

• The measured ratios Tsiin/rsiiv provide an upper 
limit T < 10 4 9 K on the temperature of the bulk of 
the gas responsible for the Si IV absorption. The 
measured rsiin/rsiiv versus Tsav and z are roughly 
in accord with, but ~ 0.15 dex lower than, the pre- 
dictions of the simulations. This may be an indi- 
cation that a small fraction of the observed gas is 
at higher temperature than in the simulations. 

• Our inferred [Si/C] is ~ 0.2 dex higher than that 
predicted by Population II, Type II supernova 
yield calculations based on a standard IMF up to 
~ 40 M©, and that observed in metal-poor stars, 
and much higher than predicted for the yields 
of massive (M < 40 M©) Population III stars. 
High [Si/C] values can, however, be obtained from 
an IMF that includes supermassive Population III 
stars exploding as pair-instability supernovae. Al- 
ternatively, we could conclude that our fiducial 
UVB is too hard; however, a UVB model signifi- 
cantly softer than our fiducial one leads to unphysi- 
cal results such as negatively evolving C metallicity 
(see Paper II). 

• Combining our measured [Si/C] with the measure- 
ments of [C/H] of paper II, we find that the Lya 
forest contributed [Si/H]= —2.0 to the global sili- 
con abundance for the QG UVB, or f^si — 3.2 x 
10" 7 . This would constitute « 30 - 50% of the 
ex pected S i prod uction by z = 2.5 as estimated 
bv lPettinil 



We are grateful to the ESO Archive for their effi- 
cient work. Without their help this work would not 
have been possible. AA and JS gratefully acknowledge 
support from the W. M. Keck foundation, and JS ac- 
knowledges support from NSF grant PHY-0070928. TT 
thanks PPARC for the award of an Advanced Fellow- 
ship. WLWS acknowledges support from NSF Grant 
AST-0206067. MR is grateful to the NSF for grant 
AST-00-98492 and to NASA for grant AR 90213.01-A. 
We thank Jason Prochaska and Stan Woosley for help- 
ful conversations. This work has been conducted with 
partial support by the Research Training Network "The 
Physics of the Intergalactic Medium" set up by the Eu- 
ropean Community under the contract HPRN-CT2000- 
00126 RG29185 and by ASI through contract ARS-98- 
226. Research has been conducted in cooperation with 

25 Note that because we can only obtain percentiles in 
T "Silv/ T ClV near and above the median, our data only constrains 
the upper half of the distribution. 



14 



Aguirre et al. 



Silicon Graphics/Cray Research utilising the Origin 2000 supercomputer at DAMTP, Cambridge. 

REFERENCES 



Aguirre, A., Schaye, J., & Thcuns, T. 2002, ApJ, 576, 1 (Paper I) 
Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 
197 

Aracil, B., P etitjean, P., Pi chon, C, & Bergeron, J. 2003, ApJ, 

submitted; astro-ph/0307506 
Boksenberg, A., Sargent, W.L.W., & Rauch, M. 2003, ApJS, 

submitted; astro-ph/0307557 
Carretta, E., Gratton, R. G., & Sneden, C. 2000, A&A, 356, 238 
Chieffi, A. & Limongi, M. 2002, ApJ, 577, 281 
Cowie, L. L. & Songaila, A. 1998, Nature, 394, 44 
Cowie, L. L., Songaila, A., Kim, T., & Hu, E. M. 1995, AJ, 109, 

1522 

Croft, R. A. C, Weinberg, D. H., Katz, N., & Hernquist, L. 1997, 
ApJ, 488, 532 

Dave, R., Hellsten, U., Hernquist, L., Katz, N., & Weinberg, D. H. 

1998, ApJ, 509, 
Depagne, E. et al. 2002, A&A, 390, 187 

D'Odorico, S., Cristiani, S., Dekker, H., Hill, V., Kaufer, A., Kim, 

T., & Primas, F. 2000, Proc. SPIE, 4005, 121 
Ellison, S. L., Lewis, G. F., Pettini, M., Chaffee, F. H., & Irwin, 

M. J. 1999, ApJ, 520, 456 
Ellison, S. L., Songaila, A., Schaye, J., & Pettini, M. 2000, AJ, 120, 

1175 

Ferland, G. J. 2000, Revista Mexicana de Astronomia y Astrofisica 

Conference Series, 9, 153 
Ferland, G. J., Korista, K. T., Verner, D. A., Ferguson, J. W., 

Kingdon, J. B., & Verner, E. M. 1998, PASP, 110, 761 
Haardt, F. & Madau, P. 2001, to be published in the proceedings 

of XXXVI Rencontres de Moriond, astro-ph/0106018 
Heger, A. & Woosley, S. E. 2002, ApJ, 567, 532 



Kim, T.-S., Cristiani, S., & D'Odorico, S. 2002, A&A, 383, 747 
McWilliam, A., Preston, G. W., Sneden, C, & Searle, L. 1995, AJ, 
109, 2757 

Norris, J. E., Ryan, S. G., & Beers, T. C. 2001, ApJ, 561, 1034 
Pieri, M., & Haehnelt, M. 2003, MNRAS, submitted; 

|astro-ph/0308003 | 

Pettmi, M. 2003, astro-ph/0303272 

Ryan, S. G., Norris, J. E., & Beers, T. C. 1996, ApJ, 471, 254 
Savin, D. W. 2000, ApJ, 533, 106 
Schaye, J. 2001, ApJ, 559, 507 

Schaye, J., Aguirre, A., Kim, T., Theuns, T., Rauch, M., & Sargent, 

W.L.W. 2003, ApJ, 596, 768 
Schaye, J., Rauch, M., Sargent, W. L. W., & Kim, T. 2000a, ApJ, 

541, LI 

Schaye, J., Theuns, T., Leonard, A., & Efstathiou, G. 1999, 

MNRAS, 310, 57 
Schaye, J., Theuns, T., Rauch, M., Efstathiou, G., & Sargent, 

W. L. W. 2000b, MNRAS, 318, 817 
Schippers, S., Miiller, A., Gwinner, G., Linkcmann, J., Saghiri, 

A. A., & Wolf, A. 2001, ApJ, 555, 1027 
Songaila, A. 1998, AJ, 115, 2184 
Songaila, A. 2001, ApJ, 561, L153 
Songaila, A. & Cowie, L. L. 1996, AJ, 112, 335 
Theuns, T., Viel, M., Kay, S., Schaye, J., Carswell, R. F., & 

Tzanavaris, P. 2002, ApJ, 578, L5 
Umeda, H. & Nomoto, K. 2002, ApJ, 565, 385 
Vogt, S. S. et al. 1994, Proc. SPIE, 2198, 362 
Woosley, S. E. & Weaver, T. A. 1995, ApJS, 101, 181 



